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Abstract 

We treat the stability of determining the boundary impedance of an ob- 
' stacle by scattering data, with a single incident field. A previous result 

! by Sincich (SIAM J. Math. Anal. 38, (2006), 434-451) showed a log sta- 

bility when the boundary of the obstacle is assumed to be C 1,1 -smooth. 
We prove that, when the obstacle boundary is merely Lipschitz, a log- log 
type stability still holds. 
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1 Introduction 

In this note we consider the stability of the determination of an impedance 
coefficient on an inaccessible boundary from scattering data arising from a single 
incident wave. 
More precisely we consider 



( Au + k 2 u = 0, inR 3 \U, 

1 &i+i\u = o. (man. 

o 

where u = u s + exp(ifcx • u), u) G § 2 and the scattered field u s satisfies the usual 
Sommerfeld radiation condition 

( du s \ 

■ lim r — iku s = 0, r = \x\ . (1.2) 

r— >oo V OT } 

The scattered field has a well-known asymptotic behavior at co 

«•(*) = 5^*0 | Uoo(J)+0 (i)} (1.3) 

as r tends to oo, uniformly with respect to x = -r^j and where Uqo is the so-called 
far field pattern of the scattered wave (see [TO])- 
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We deal with the inverse scattering problem of determining A = X(x) when the 
far field Uqo is known for a given incident direction uj. 

The stability for this problem was treated in [TO] , and we refer to the bibliogra- 
phy therein for further references. 

The result in [TO] can be summarized as follows. If the unknown impedance A 
is Lipschitz and the boundary of the scatterer D is C 1 ' 1 , then A depends on Uoo 
with a modulus of continuity of logarithmic type (with a single log !). 

The purpose of this note is to investigate the stability when the regularity 
assumption on the scatterer is relaxed to Lipschitz. The requirement of the 
C l > regularity at the boundary was mainly due to the fact that such regularity 
is needed in order to obtain a doubling inequality at the boundary for \u\ 2 [TO1 
Theorem 4.6]. We refer to [T| for a thorough discussion of this topic. In brief, 
such a doubling inequality is needed because the impedance can be obtained 
from the total field u (which is uniquely, and stably, determined by the far field 
pattern) by computing 

A * ^ U 
u dv 

Consequently, it is necessary, in order to estimate A, to have a control on the 
rate of vanishing of u on dD. The doubling inequality implies that such a rate 
is at most algebraic in the L 2 -average sense, that is 




for some constant K > 0. Here A r (xo) = B r (xo) H dD and B r (xo) is the ball 
of radius r centered at xq £ dD. From this bound on the vanishing rate of |it| 2 
the stability for A with a single log follows. Note that in fact, in [TO] a slightly 
different route is followed, which involves the notion of Muckenhoupt weights. 
It is an open problem if the doubling inequality at the boundary holds true 
when we relax the regularity assumption on dD (see again fj]). Nevertheless, it 
is expected that the order of vanishing of u\qd can be still controlled somehow. 
In fact we shall prove that, assuming dD Lipschitz, u\q D has a rate of vanishing 
controlled in an exponential fashion as follows 

/ \uf>e-«r- K . (1.4) 

JA r (x a ) 

With the aid of such estimate we are able to prove a stability result of a weaker 
type, namely with a log-log type modulus of continuity. 

We remark that this phenomenon confirms a scheme that was already known for 
different, but related inverse boundary value problems [2]. In that paper, inverse 
problems with unknown boundary for elliptic equations were treated and, on the 
unknown boundary a condition of Dirichlet or Neumann type was prescribed. 
It was shown that, in case that the unknown boundary is smooth enough to 
have a doubling inequality at the boundary, then the stability was of log type. 
Instead with weaker regularity assumptions only a log-log type estimate could 
be obtained. See for instance the proof of [U Theorem 2.1]. This scheme has 
also been confirmed with different inverse boundary problems |14[ [5] . 
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Let us mention also that, recently, a similar task to the present one was at- 
tempted in the preprint [7J by Bellassoued, Choulli and Jbalia. The main result 
in that paper also consists of a log-log stability, but of local type only. That 
is, the stability is proved only when the error on the unknown impedance is 
a-priori known to be small. Furthermore the stated regularity assumptions on 
the boundary have some unclear aspects. 

The plan of the paper is as follows. In Section[5]we introduce notation, the main 
assumptions and state the main stability result, Thcorcm l2.ll Section [3] collects 
more or less well-known estimates on the direct problem , wc conclude it by 
stating a rather more delicate result, Theorem 13.41 in which the boundary trace 
of the field u is estimated in H , The proof is contained in the following Section 
2J It relies on the well-known Rellich's identity, the main theorem of this section 
being Theorem 14. II Section [5] explores the stability of the determination of the 
scattered field in terms of the far field data. We review some known results, 
and we present some new ones, which rely on the analysis developed in the 
previous Section [H see in particular Corollaries 15.31 and 15.41 where the stability 
for negative norms of the normal derivative of u is obtained. The occurrence 
of negative norms appears to be necessary here because, being the boundary 
only Lipschitz, the normal derivative of u may be nonsmooth. In Section [5] 
we obtain, Theorem 16.11 the local lower bound on the vanishing rate of u on 
dD announced above in (jl.4|) . This is achieved though a quantitative estimate 
of unique continuation which is by now well-known as a Lipschitz Estimate of 
Propagation of Smallness, see for instance [5J[33]. In Section [7J we present an 
interpolation inequality between a weighted L^-norm and a Holder norm, when 
the weight satisfies a bound on its vanishing rate of the type (|1.4[) , Proposition 
17.11 We finish the proof of Theorem 12.11 in the final Section [8] 

2 Main assumptions and results 

2.1 Main hypotheses and notation 
Assumptions on the domain. 

We shall assume throughout that D is a bounded domain in R 3 , that is, for a 
given d > 0, we require diamD < d. Also we require that R 3 \ D is a connected 
set, and that the boundary dD is Lipschitz with constants ro, M. More precisely, 
for every xq G dD, there exists a rigid transformation of coordinates under 
which, 



where x G R 3 



x = 



Df}B ro (x ) = {(x',x 3 ) : x 3 > 7 (x')} 
= (x', x 3 ), with x' G R 2 , x 3 G R and 



(2.1) 



7 : B' ro {xo) Cl 2 4l 



satisfying 7(0) = and 



H7||c<m(b;. Oo)) < Mr o, 



where we denote by 



|7||z»(b; o Oo)) + r o su P 

x,y£B' ro (z ) 



7p) - 7(2/) 

\x - y\ 
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and B'(xq) denotes a ball in K 2 . 

For the sake of simplicity we shall assume that G D. 

Fixed R> d, p G (0,r ) and xq G dD, let us define the following sets 



D + =W\D, (2.2) 

D+=B R (Q)nD+, (2.3) 

E p = {x e R n \ < diat(x,D) < p)}, (2.4) 

D p = D U E p , (2.5) 

T p (x ) = B p (x )\D, (2.6) 

A p (x ) = B p (x )ndD. (2.7) 

A priori information on the impedance term. 

Given Ao > 0, we assume that the impedance coefficient A belongs to C 0,1 (dD, R) 
and is such that 

\{x) > A > (2.8) 

for every x G dD. Moreover we assume that, for a given constant A > 0, we 
have that 

l|A||co.i(8D) < A- (2-9) 

From now on wc shall refer to the a priori data as to the following set of 
quantities, d, ro, M, Ao, A, k. 

In the sequel we shall denote with rj(t) a positive increasing concave function 
defined on (0, +oo), that satisfies 

V(t) < C(log(t))-<>, for every < t < 1 , (2.10) 

where C > 0, > are constants depending on the a priori data only. 

2.2 The main result 

Theorem 2.1 (Stability for A). Let X u A 2 satisfy (|2~%]l . ([23| . Let it,, i = 1,2, 
6e t/ie weafc solutions to the problem (|1.1[) wit/i A = A-; respectively and let Ui.oo 
fee t/iezr respective far field patterns. If for some e > 0, we have 

||«l,oo - U2,oo||i»(8Bi(0)) < £j (2- 11 ) 

t/ien 

l|Ai-A 2 || L » (a£)) < i? 077(e), (2.12) 
where r\ is given by (|2.10p . 

3 The direct problem 

Let us introduce the following space 

Hl QC {D+) = {»£ V{D+) : v\ D + G H\D+), for every i? > s.t. 2} C B^(0)} 
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where V is the space of distribution on D + . 

A weak solution to the problem is a function u = exp (ikiu ■ x) + u s , where 
u s G H^ qc (D + ) is a weak solution to the problem 

Au s + k 2 u s = 0, inD+, 
du s 

— h i\(x)u s = ~~q~ exp (ikto ■ x) — i\{x) exp (ikuj ■ x), on 3-D, /g 

V ( du s V „ \ 

lim r _ ) . 00 r I (rx) — iku s (rx) j = 0, uniformly in i. 

where is the inward unit normal to D. 

Lemma 3.1 (Well-posedness). The problem (|3.ip has one and only one weak 
solution u s . Moreover, for every R> d, there exists a constant Cr > depend- 
ing on the a priori data and on R only, such that the following bound holds 

\\u s \\ h1{d+r) <C r . (3.2) 

Proof. For the proof we refer to Lemma 3.1 in [TO] , see also [5] for a previous 
related result. □ 



Theorem 3.2 (C a regularity at the boundary). Let u be the weak solution 
to (jl.ll) . then there exists a constant a, < a < 1, such that for every R > d 

u G C a (-DftJ ■ Moreover, there exists a constant Cr > depending on the a 
priori data, on R only, such that 

\\u\\ ca(5f) <C R . (3.3) 

Proof. This is a more or less standard regularity estimate up to the boundary. 
The Moser iteration technique fits to this task. Details can be found in [TSJ 
Lemma 3.3]. Note also that the arguments used there only require the Lipschitz 
regularity of dD. □ 



Corollary 3.3 (Lower bound). There exists a radius Rq > depending on 
the a priori data only such that 

\u(x)\ > - for any \x\ > Rq . (3.4) 

Proof. The proof relies on the same arguments discussed in [191 Corollary 3.3]. 
□ 



Theorem 3.4 (^(dD)). Let u be as in Theorem\K^ then 

\\u\\ H HdD) < C . (3.5) 
The proof shall be given in the next section. 
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4 Estimates at the boundary. 

We begin by recalling the well-known fact that there exists po G (0, ro] depend- 
ing only on M,ro,d and k only such that, for every p G (0, po] the following 
coerciveness condition is satisfied 




for all v G Hq (E p ), see for instance [TTJ Lemma 8.4]. In what follows we shall 
fix p G (0, pol- 



Theorem 4.1. Let v G H 1 (E p ) be a solution to 

Av + k 2 v = in E p . 
If its trace v^gjj G H 1 (dD) then ^: gD G L 2 (dD) and we have 



dv 
dv 



< C 



L 2 (dD) 



{\\W T v\\ 



2 

L 2 (dD) 



Hi 



(4.2) 



(4.3) 



Conversely, if §f;, dD G L 2 (dD) then v\o D G H 1 (dD) and 



It 2 (a 





dv 




dv 



L 2 (dD) 



(£'„ 



(4.4) 



Here Vtv denotes the tangential gradient ofv on dD and C depends on M, ro, d, p 
and k only. 

Proof. A priori inequalities of this sort were first proven by Payne and Wein- 
berger [15| , [16] . A proof when k = and with Lipschitz boundary is due to 
Jcrison and Kenig [13] . The underlying tool for such estimates relies on the 
celebrated Rellich's identity [T7]. The adaptation to the case of the Hclmholtz 
equation (k ^ 0) can be obtained as follows. 

Recalling the local graph representation of dD introduced in (]2.ip . let us define, 
for every r G (0, 2po], the half cylinder 



C+ = {x = (x',x 3 ) G M 3 ||a;'| < r, 1 (x 1 ) < x 3 < Mr} 
and let us denote 

A r = {x = (x',x 3 ) G M 3 ||a;'| < r, x 3 = J (x')} . 
First we prove local estimates of the following form 



dv 
dv 



< C 



|vh + / Iv^r + K 



(4.5) 
(4.6) 

(4.7) 




G 



for every r E (0,Po\ an d C depends on M,ro,r and k only. These estimates 
are proven in the case k = in 4, Proposition 5.1]. The extension to the 
case k ^ can be either obtained by a straightforward adaptation of the argu- 
ment in [H Section 5], or else by introducing an auxiliary variable t and posing 
V(x,t) = v(x)e kt which is harmonic in C^ r x K. Writing down the analogs in 
higher dimension of (|4.7[) and (|4.8I) for V, we readily deduce the above stated 
inequalities for v. The final step consists of covering dD with neighbourhoods 
A n with r small enough so that the corresponding half cylinders C^ r remain 
within Ei p . The number of the neighbourhoods A r needed to cover dD can be 
estimated in terms of M, ro, d, p, see [H Proof of Prop. 3.3]. □ 

Proof of Theorem 13.41 Immediate consequence of (|4.4[) in combination with 
the impedance condition in (jl.l[) and Lemma l3.ll □ 



We shall denote by $ E C°° (M. n ) a mollification of the characteristic function 
XDp/2 such that O<0<M = lin D p/3 , i? = in M"\£ 2p/3 , \D m ti\ < C m p~ m 
in R™, where C m depends on m, d and ro. 

Lemma 4.2. Let u E H l (E p ) n C° (i? P ), &e a solution to 

Au + k 2 u = in E n . 



We have 



du 
dv 



(4.9) 
(4.10) 



<C\\u\\ L ^ (Ep) 

H- 1 (dD) 

where C depends on M, rg,d, p and k only. 

Proof. In view of (|4.1[) . for any C E H 1 (dD) we can consider the unique 
solution tp E H 1 (E p ) to the Dirichlet problem 



Aip + k 2 ip = in E p , 

if = ( on dD , 
tp = on <9£ p \ dD . 



Moreover we have 



(4.11) 



(4.12) 



IMIffi(.E p ) - C|lCllffi/2(a.D) 

with C > only depending on M, r$,d, p and fc. Let § the previously introduced 
cutoff function and denote ip = "dtp. The Green's identity gives 



, du f dip 

dD OV Jg D dv 



u (Aip + k 2 ip) 



hence 



/ < / + f u (2Vi9 • Vip + (Atf + k 2 -d) ip) 

JdD 9v J dD dv J Ep 

applying (|4.3p to <p and taking into account (|4.12j) we get 



(4.13) 



(4.14) 



- du I 



< 



I JdD C du 

< C (\\u\\ Lao(Ep) K\\hH8D) + \Hl°°(e p ) \\<P\\b^B p ) 



< 



<C\\u\ 



L°°(E P 



IKIIhi 



(4.15) 
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and the thesis follows by duality. 



□ 



5 Stability for the scattered field 

Lemma 5.1 (From the far field to the near field). Let Uj,iij i00 , i = 1,2, 
be as in Theorem \2.1l Suppose that, for some e, < e < 1, (|2.11j) holds, then 
there exist a radius i?i > and a constant C > 0, depending on the a priori 
data only, such that 

||ui-M2|U2(b Ri+1 (o)Vb Ri (o)) < Ce a{£ \ (5.1) 
where a(e) is defined as 

a(t) = : — t, — ~7 — 77 r . (5.2) 

Proof. For the proof we refer to Lemma 4.1 in [19], which is based on the sta- 
bility results for the near field achieved by Isakov in [T2] and further developed 
by Bushuyev in [5]. □ 



Theorem 5.2 (Stability at the boundary). Let Ui,Ui j00 , i = 1,2, be as in 
Theorem \2.1[ We have that, if for some e > 0, (|2.11[) holds, then 

\\ui - u 2 \\ L ^(dD) < ??(e) , (5.3) 

where r\ is given by (|2.10[) . with a constant C > depending on the a priori 
data only. 

Proof. The arguments in need few adjustments. The main additional tool 
is a global estimate of propagation of smallness [SJ Theorem 5.3] which enables 
to achieve 

- u 2 \\ L ^ D + ) < r){s) . (5.4) 

The C a bound obtained in Theorem 13.21 allows, by interpolation, to conclude. 

□ 



Corollary 5.3 (H 1 bound). Under the same hypothesis of Theorem \3.2\ we 

have that 



dui 
dv 



du 2 
dv 



H- l (dD) 

Proof. This is an immediate consequence of (|5.3[) and of Lemma l4~2l 



(5.5) 
□ 
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Corollary 5.4. Under the same hypothesis of Theorem \3.2\ we have that 



dui 


du 2 




dv 


dv 


H 



<ri{e) 



(5.6) 



2 (dD) 



Proof. The result follows by interpolation with the aid of the impedance con- 
dition in (fTTj) and of Theorem [321 □ 



6 The estimate of Lipschitz propagation 
of smallness 

Let us denote r\ = min {ro, ^}. 

Theorem 6.1. Let u be the weak solution to flLTJ) . For every r, < r < r\, 
and for every xq £ dD we have that 



( \u\ 2 > exp(-Kr- K ) 



IA r {x ) 

where K > only depends on the a priori data. 
The proof shall stem from the two Lemmas below. 

Lemma 6.2. Let u be the weak solution to For every r, < r < ro and 

for every xq £ dD we have that 

(7 W 2 ) >C [ M 2 (6.2) 

\Jaaxo) J Jr,(x ) 

where C > 0, 0</3<^ are constants depending on the a priori data only. 

Proof. By the arguments in [5j Theorem 1.7] concerning the well-known local 
estimate for the Cauchy problem and by the a priori bounds achieved in Theorem 
and Theorem 13.41 we get that for any < r < ro and for any xo £ dD 

IMU^rv^o)) < 

* C (lHU(A r( * 0) ) + IMh- Wo)))' OMMr^o))) 1 - 5 (6-3) 



where C>0,0<<5<1 arc constants depending on the a priori data only. Note 

,1-5 

l£ 2 (rv(z )) 



that the last factor 1 1 u 1 1 ^ 2 ^ , is bounded by a constant in view of Theorem 



We recall that for any < r < ro and for any xo £ dD the following interpolation 
inequality holds 

ll U ll^(A,(, ))^ C ll"lli(A,.(, ))ll-lli(A r( , 0) ) (6-4) 

where C > depends on the a priori data only. 

Finally by the impedance condition on dD, by the above interpolation inequal- 
ity and by the a priori bound in Theorem l3.4l the theorem follows with f3 = |. □ 
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Lemma 6.3. Let u be as above. For every r, < r < ^ and for every xq G dD 
we have that 



\u\ 2 > Cexp(-fcir- fc2 ) (6.5) 

r r (x ) 

where C, fci , &2 > are constants depending on the a priori data only. 

Proof. For the proof we mainly refer to [Ml Proposition 3.1], where the authors 
achieved a Lipschitz propagation of smallness result for the L 2 norm of Jacobian 
matrix of the solution to the Lame system. Nevertheless, the same arguments 
go through for the L 2 norm of the solution to the Hclmholtz equation as well. 
We provide below a sketch of the proof. 

Let x £ r r (:co) be such that Bn{x) C Y r (xo). Now, by adapting the techniques 
developed in |14l Proposition 3.1], which are mostly based on a standard prop- 
agation of smallness and the iterated use of the three spheres inequality along a 
chain of balls centered on the axis of a cone and with increasing radii (see also 
[3J Section E.3]), we find that 

t| 2 > Ccxp(-k ir ~ k2 ) I \u\ 2 (6.6) 



Br(x) J D+ 



-211,, 



where Ro > is the radius introduced in Corollary 13.31 
By combining the lower bound stated in (|3.4[) and the obvious inequality 
L , \ lul 2 > f D \u\ 2 we obtain the desired estimate. □ 



Proof of Theorem 16.11 From the above Lemmas we deduce 

\u\ 2 >ctexp(-^r- fc2 ) . (6.7) 



It is an elementary fact that we may find K > 0, only depending on C, ki, &2, (3 
such that (|6.1|) follows. □ 



7 A weighted interpolation inequality 

Proposition 7.1. Given M, K > 0, let w > be a measurable function on dD 
satisfying the conditions 

\\w\\L°°(dD) < M (7.1) 

and 

ll w; l|L 2 (A r (a:)) > ex p(~ Kr~ K ) for every x S dD and r 6 (0,ri). (7.2) 
Let f e C a {dD) be such that 

\f{x)~ f(y)\<E\x-y\ a for every x, ye 3D. (7.3) 

If 

I \f\w<e 

JdD 
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(7.4) 



the 



ll/l 



< Er) 



\ L °°(dD) _ -., 

where 77 satisfies (|2.10[) with constants only depending on M, K,ro,a. 
Proof. By (|7.1[) we have trivially 



(7.5) 



/ w > M" 1 exp(~2Kr~ K ) for every x E 3D and r E (0, n). (7.6) 

JA r (i) 

Now fix £ e <9D such that |/(£)| = For every r > and ir G A r (£) 

we have 

1/(01 < \f(x)\+Er a , (7.7) 
by multiplying both sides of (|7.7|l by w and integrating over A r (£) we get 



1/(01 / ™< / H/l+^ a / 



(7.8) 



hence 



II/I|l«(S£») 



< 



f*r® W 

< Mexp(2Kr~ K )e + Er a . 



(7.9) 



Now, if e < E and r = {AK) 1 / K \ log %\ X I K < n, then we fix this value of r in 
(17.91) and the thesis follows. The remaining cases are trivial. □ 



8 Conclusion 

Proof of Theorem I2TT1 

By a standard interpolation estimate we deduce that 



||ui(Ai - A 2 )|| L 2(aD) ^ C||wi(Ai - A 2 



Iff 1 (dD) 



||«i(Ai 



A2)II ff -i ( a.) 



where C > only depends on the a priori data. 
It is an elementary observation, that we have 



||ui(Ai - A 2 )||ffi(aD) < II Ai - A 2 ||c o ' 1 (aD)ll w i||H 1 (aD) < C 

where C > is a constant depending on the a priori data only. 
On the other hand, by the impedance condition 



duj 
dv 



+ iXiUi = on dD 



(8.2) 



(8.3) 



for i = 1 , 2 we also have that 



ki(Ai-A 2 )|| 1 < 





du 2 




dv 


dv 


H 



C\\U! 



H~2(dD) 



u 2 \\ 1 . (8.4) 
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By combining the results in Theorem 15.21 and in Corollarv l5.4l we get 

\\u 1 (\i-\ 2 )\\ H - i(gD) <r?( £ ) . (8.5) 

where C > is a constant depending on the a priori data only. 
Hence, by the inequality (|8.1I) we arrive at 

||«i(Ai-A 2 )|| £a (8D) <»/(e) ■ (8.6) 

Now, we apply Proposition 17.11 with w = \ui\ 2 and / = (Ai — A2) 2 and the 
conclusion follows. □ 
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